The pair annihilation reaction D + D — > in disordered media and conformal invariance 
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The raise and peel model describes the stochastic model of a fluctuating interface separating a 
substrate covered with clusters of matter of different sizes, and a rarefied gas of tiles. The stationary 
state is obtained when adsorption compensates the desorption of tiles. This model is generalized 
to an interface with defects (D). The defects are either adjacent or separated by a cluster. If a 
tile hits the end of a cluster with a defect nearby, the defect hops at the other end of the cluster 
changing its shape. If a tile hits two adjacent defects, the defect annihilate and are replaced by a 
small cluster. There are no defects in the stationary state. This model can be seen as describing 
the reaction D + D — > 0, in which the particles (defects) D hop at long distances changing the 
medium and annihilate. Between the hops the medium also changes (tiles hit clusters changing 
their shapes). Several properties of this model are presented and some exact results are obtained 
using the connection of our model with a conformal invariant quantum chain. 

PACS numbers: 05.50.+q, 47.27.eb, 05.70.-a 



I. INTRODUCTION 

The raise and peel model 0, [|J which is a stochas- 
tic model of a fluctuating interface is, to our knowledge, 
the first example of a stochastic model which has the 
space-time symmetry of conformal invariance. This im- 
plies that the dynamic critical exponent z = 1 and cer- 
tain scaling properties of various correlation functions are 
known. This model was extended in order to take into 
account sources at the boundaries @, 1, H keeping con- 
formal invariance. In all these cases, the stationary states 
have magic combinatorial properties. 

In the present paper we describe another extension of 
the raise and peel model keeping conformal invariance 
(see also Appendix A and B) by introducing defects on 
the interface. These defects hop at long distances in 
a medium which is changed by the hops. Between the 
hops the medium also changes. Finally when two defects 
touch, they can annihilate. The stationary state is the 
one as the original raise and peel model with no defects. 

The whole process can be seen as a reaction D+D — » 0, 
where D is a defect, taking place in a disordered un- 
qucnchcd medium. 

In Sec. II we describe the model. Like the raise and 
peel model [![, the present model comes from consider- 
ing the action of a Hamiltonian expressed in terms of 
Temperley-Lieb generators on a vector space which is a 
left ideal of the Temperley-Lieb algebra. The ideal can 
be mapped on graphs which constitute the configuration 
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space of the model. We shortly review in Appendix A 
the mathematical background of the model and refer for 
details to Refs. 0,1]. 

In Sec. Ill, using Monte Carlo simulations, we de- 
scribe the long range hopping of defects and give the 
Levy flights probability distribution. 

In Sec. IV, again using Monte Carlo simulations, start- 
ing with a configuration which consists only of defects, 
we study the variation in time t of their density for a 
lattice of size L. We obtain the scaling function which 
gives the number of defects in terms oit/L and show how 
conformal invariance gives some of its properties. In the 
thermodynamic limit the density decreases in time like 
1/t as is expected since in a conformal invariant theory 
time and space are on equal footing. In Sec. V we present 
our conclusions. 



II. THE RAISE AND PEEL MODEL WITH 
DEFECTS. 

We consider an interface of a one-dimensional lat- 
tice with L + 1 sites. An interface is formed by at- 
taching at each site a non-negative integer height hi 
(i = 0,1,..., L). We take h = h L = 0. If for two 



consecutive sites j and j + 1 we have hj 



0, on 



the link connecting the two sites we put an arrow called 
defect (see Fig. 1). For the remaining sites, the heights 
obey the restrict solid-on-solid (RSOS) rules: 



±1, 



hi > 0. 



(1) 
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A domain in which the RSOS rules are obeyed {hj = 
hi = 0, hk > 0, j < k < 1} is called a cluster. There are 
3 clusters and 3 defects in Fig. 1 (L = 21). There are 
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FIG. 1: One of the configurations for L — 21 (22 sites). 
There are 3 defects (arrows on the links) and 3 clusters. Also 
shown are 5 tiles (tilted squares) (a)-(e) belonging to the gas. 
When a tile hits the surface the effect is different in the five 
cases. 

([I//2]) P 0SS1 ble configurations of the interface (we denote 
by [x] the integer part of a;) . 

There is a simple bijection between the configurations 
of interfaces and defects, where Fig. 1 is an example, 
and ballot paths |4|]. A ballot path is obtained if one 
follow the RSOS rules ([T]), take h L = but let h free 
(0 < h,L < L) . This fact was used in Q to define another 
stochastic model than the one described below. In the 
case L = 4, the six possible configurations are shown in 
Fig. 2. The configuration shown in Fig. 2b has 2 defects 
and one cluster, while there are no defects in Fig. 2f. 

, t , t , t , t . /\ t , t , „ 



tiles (this is an odd number). Similarly, if the 
slope is negative (/li+i < hi < the tile is 

reflected after triggering the desorption of a layer 
of tiles belonging to the segment (hj > hi = hi-i, 
j + i-b + l,...,i-l). 

d) The tile hits the right end of a cluster hj > h(i — 
c) = h(i) = (j = i-c+l,...,i-l)andh(i+l) = 0. 
The link (i, i + 1) contains a defect ("d" in Fig. 1). 
The defect hops on the link (c, c+1) after triggering 
the desorption of a layer of tiles (hj — ► hj — 2, j = 
i — c+ 1 , . .. , i — 1) and the tile is adsorbed producing 
a new small cluster = h i+ i = 0, hi = 1) (see 

Fig. 3). If the defect is at the left end of a cluster, 
the rules are similar, the defect hops to the right 
after peeling the cluster, and a new small cluster 
appears at the end of the old one. 
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FIG. 3: The new profile after the tile "d" in Fig. 1 has hit 
the surface at the right end of a cluster. The defect hops to 
the left end of the cluster, peeling one layer and a new small 
cluster appears at the right of the old cluster. 

e) The tile hits a site between two defects = 
hi = h i+ i = 0). This is the case "e" in Fig. 1. The 
two defects annihilate and in their place appears a 
small cluster (hi-i = hi+i = 0,hi = 1). See Fig. 4. 




FIG. 2: The six configurations for L = 4 (5 sites). In the 
stationary state only the RSOS configurations e) and f) occur. 

We consider the interface separating a film of tiles 
(clusters with defects) from a gas of tiles (tilted squares) . 
The evolution of the system (Monte Carlo steps) is given 
by the following rules. With a probability Pi — a 
tile from the gas hits site i (i = 1, . . . , L— 1). As a result 
of this hit, the following effects can take place: 

a) The tile hits a local maximum of a cluster (" a" in 
Fig. 1). The tile is reflected. 

b) The tile hits a local minimum of a cluster ("b" in 
Fig. 1). The tile is adsorbed. 

c) The tile hits a cluster and the slope is positive 
(h l+1 > hi > fti_ x ) ("c" in Fig. 1). The tile 
is reflected after triggering the desorption of a 
layer of tiles from the segment (hj > hi = 

j = i + 1, . . . , i + b — 1), i.e., hj — > hj — 2, 
j = i + 1, + b — 1. The layer contains 6 — 1 
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FIG. 4: The new profile after the tile "e" has hit the surface 
between two defects. The defects have disappeared and in 
their place one gets a new small cluster. 

To sum up. The defects (D) hop non-locally in a dis- 
ordered (not quenched) medium which changes between 
successive hops (local adsorption and nonlocal desorption 
take place in the clusters). During the hop the defect 
peels the cluster and therefore also changes the medium. 
The annihilation reaction D + D — > is local. If one 
starts the stochastic process with a certain configuration 
(for example, only defects as in Fig. [Sk), due to the an- 
nihilation process, for L even all the defects disappear 
and in the stationary state one has only clusters (RSOS 
configurations). The properties of the stationary states 
have been studied elsewhere [H, 0- In the case L odd, 
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in the stationary states one has one defect. In the next 
section we are going to see how this defect hops and will 
observe that the defect behaves like a random walker per- 
forming Levy flights. This will help us understand the 
annihilation process D + D — ► described in Sec. IV. 
The rules described above were obtained by using a rep- 
resentation of the Temperley-Lieb algebra in a certain 
ideal 0, H d(see Appendix A). The finite -size scaling 
of the Hamiltonian eigenspectrum is known from confor- 
mal field theory (see Appendix B), therefore the physical 
properties of the model can be traced back to conformal 
invariance. 



III. 



THE RANDOM WALK OF A DEFECT 



Before discussing the annihilation reaction of defects, 
it is useful to understand how defects hop. The simplest 
way to study the behavior of defects is to take the sta- 
tionary states in the case L odd when we have only a sin- 
gle defect. Although there is a lot of information about 
these stationary states coming from combinatorics 0, 
and Monte Carlo simulations p|, the results we present 
here are new. 

One asks what is the probability P(s) for a defect to 
hop, at one Monte-Carlo step, at a distance s (we assume 
L very large). We first see if on physical grounds, one 
can't guess the result. Let us assume that the defect 
behaves like a randon walker and that P(s) describes 
Levy flights @, i, H, [nj. This implies that for large 
values of s we have: 



P(s) 



1 



ll + cr ■ 



(2) 



If the random walker starts at a point x — (for example 
in the middle of the lattice), at large values of t, the 
dispersion is [l(J: 



(3) 



In a conformal invariant model, one has no other scales 
but the size of the system, space and time are on equal 
footing therefore one has to have <r = 1. 

In Fig. 5 we show P(s) as obtained from Monte Carlo 
simulations for systems of different sizes. One notices a 
data collapse for a large domain of s. A fit to the data 
for the largest lattice (L = 4095) gives, for large s: 



P(s) 



2.25 

I „ 12.06 : 



(4) 



in agreement with what we expected. 



IV. THE DENSITY OF DEFECTS AT LARGE 
TIMES 

We are now going to study the number of defects 
Nd(t, L) as a function of time and lattice size taking at 




ln(lsl) 

FIG. 5: (Color online) The probability P(s) for a defect to 
hop at a distance s in units of lattice spacing. Monte Carlo 
simulations were done on systems of different sizes. 



t = the configuration where the lattice is covered by 
defects only (like in Fig. 2a). An interesting and novel 
aspect of this study is the role of conformal invariance. 
Since there are no other scales in the system except L, 
we expect for large values of t and L: 



N d (t,L) = f(j.). 



(5) 



In Fig. 6, we show Nd(t, L) for various lattice sizes (L 
odd) . One sees a nice data collapse except for very small 
values of t/L where the convergence is slower. A similar 
(but not identical!) function is obtained for L even. 

We firstly discuss the behavior of Nd for large values 
of t/L (see Fig. 7). A fit to the data gives (L odd): 



N d [j) 



4°v 



where 



Co) 
l 

(°) 



6.75, 
8.21, 



A 



(°L 



" A o T 



(o) 



: 17.27, 
26.48. 



(6) 



(7) 



We can now compare the data obtained from the fit with 
the finite-size scaling spectrum of the Hamiltonian (see 
Appendix B, Eqs. (lEgf . (|B9| ) : 



A 



(o) 



3tt\/3 _ 
2 

3nV3 10 



.1620971 



= 27.20699 • 



(8) 



No prediction can be made about or A^' since they 
are not universal, they depend on the initial conditions. 
Notice that A { 2 o) > A{ o) as it should be since the expan- 
sion should diverge for short times where we expect 

L 
1 



do) 



N, 



(9) 



4 



L=4095 
- - L=2047 
^L=1023 
*-*L=511 




FIG. 6: (Color online) The number of defects Nd(t,L) as a 
function of t/L for several lattice sizes (L odd). At t = 0, 
Nd(t = 0,L) = L. The error bars are also shown. The fitted 
linear curve shows that the density decreases as the inverse of 
time. 
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FIG. 7: (Color online) The number of defects Nd as a function 
of t/L as in Fig. 6 zoomed on the large time domain. The 
error bars, given in Fig. [6] are not shown. 

A similar fit, done for L even (the data are shown in 
Fig. 8), gives: 



with 



2.83, A 



6.93, 



\[ e) = 2.71, 



\ { 2 e) = 16.64. 



(10) 



(11) 



We can again use the predictions of conformal invariance 
(see (jB8j) and ([B9])) and get 



3ttV3 1 
2 3 



2.72069- ••, 



A 



(e) 



37rV3 r 



16.324194- ■■ 



(12) 



to be compared with (fTl"]) 
5 




FIG. 8: (Color online) Nd as a function of t/L for large times 
for different lattice sizes (L even). The error bars are of the 
same order as in Fig. [6] 



In the small t/L domain we get for L even and odd 

0.322 



P : 



Nd 
L 



t 



(13) 



In order to find the correction term in (|13|) , we have com- 
puted Ndt 2 /L as shown in Fig. 9. We have obtained a 
straight line from which we get: 



0.322 0.334 



t 



t 2 



(14) 



This last result is the same for L even and odd. Notice 
that the correction term in (|14p is not given by the scaling 
function 

We have also computed the fluctuation of the density 
as a function of time and got: 



< p 2 > - < p > 2 _ 0.237 
< p > 2 ^ t lm ' 



(15) 



We would like to compare our results with known re- 
sults obtained for diffusion and annihilation reactions ( 
A + A -> ) with Levy flights [H OS OS 03 01- In 
one dimension, for Levy flights given by Eq. ((2]), one gets 
El: 



" for 
l£t for 



for 



a > 1 

(7 = 1 

a < 1 



(16) 



the critical dimension being d c = a . 

If one compares (TTTjl) for a = 1, as obtained in Sec. 
Ill and ([14"]) , one notices the absence of the \nt correc- 
tion. Such a term, if present, could have been seen in 
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t 



FIG. 9: (Color online) The density of defects times t 2 for 
short times. A linear fit to the data obtained for the largest 
lattice (L = 4096) gives (fl4"|. 

our simulations (one observes that for large lattices, pt 
converges to the value 0.322 from above). Logarithmic 
corrections can also appear in a conformal field theory if 
one has Jordan cells [l6| but there are no Jordan cells 
in the Hamiltonian (|B2|) given in Appendix A [13] ■ We 
believe that the discrepancy between the results of our 
model and those obtained for the reaction A + A — > 
comes from the fact that the two models have little in 
common. 



V. CONCLUSIONS 

We have presented an extension of the raise and peel 
model taking into account defects. The main property 
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APPENDIX A: THE CONNECTION BETWEEN 
TEMPERLEY-LIEB STOCHASTIC PROCESSES 
AND THE RAISE AND PEEL MODEL 

We shortly review this connection, for details see P,[l8j 
and (|. 

Consider the Temperley-Lieb semigroup algebra G de- 
fined by L — 1 generators ej (j = 1, . . . , L — 1) and the 
relations: 



e i e i±ie 
for \j 



k\ > 1, 



and the Hamiltonian 



F=^(l- ei ). 



(Al) 



(A2) 



In the basis {u> c } of the words of G (the regular repre- 
sentation of G), H is a matrix satisfying H a jj < for 
a =/= b and ^ b = 0. Such a matrix is an intensity 
matrix and defines a Markov process in continuum time 
given by the master equation 



df 



(A3) 



where P a {t) is the (unnormalized) probability to find the 
system in the state |o) at time t and the rate for the 
transition \b) — > \a) is given by —H a ^, which is non- 
negative. The Hamiltonian (|A2|) has an eigenvalue equal 
to zero. The corresponding left eigenvector (0| is trivial, 
the right eigenvector |0) gives the probabilities in the 
stationary state: 

(0\H = 0, (0| = 5>|, 

a 

H\0) = 0, |0)=VP Q |a), P Q = lim P„(t).(A4) 

* — • t— >oo 



iJ defined by (|A2|) gives a Markov process not only if it 
acts in the vector space of the regular representation but 
also if it acts in the vector space of a left ideal I because 
the generators ej map the left ideal into the left ideal: 



a = i. 



(A5) 



An easy way to define the left ideal in which we are 
interested and the action of the generators on this ideal 
is to use the language of graphs. 

The generators ej can be pictorially represented by 













1 2 


3- 


-1 3 3+1 3 


+2 L — 


1 i 



(A6) 



The elements of the ideal can be represented by links- 
defects diagrams. They can be obtained in the following 
way (see Fig. [10] for L = 4): take L sites. If a site is not 
connected to another one, draw a vertical arrow. Two 
sites can be connected by a link. The links don't cross 
each other and the arrows can't cross the links. For a 
given L the number of diagrams with m defects is: 



Cr 



L 

L-m+l - 



L 

L — m — 1 ' 



and the total number of diagrams is 



L/2 



Dd2) 



(A7) 



(A8) 



where [a;] is the integer part of x. 

Ill la / 
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1 1 



i b 




FIG. 10: The six links-defects diagrams for L — 4. The 
diagrams a-f correspond to the a-f RSOS configurations of 
Fig.H 

The action of ej on a links-defects diagram is given 
by placing the graph of ej underneath the first diagram, 
removing the closed loops and the intermediate dashed 
line. Next one contracts the links in the composite pic- 
ture. In Fig.[TT]we show the action of e2 on the diagram 
b of figure M 




1 1 



' 2 IU^ lO 

FIG. 11: The action of e-z on the diagram b of Fig. 1101 



The action of the Hamiltonian (|A2|) in the vector space 
given in Fig. fTDl is 



H 



f 3 
















-1 


2 


-1 











-1 


-1 


2 


-1 








-1 





-1 


2 











-1 





-1 


1 


-2 


Vo 











-1 


2/ 



(A9) 



Notice that H has a block triangular form. The station- 
ary state |0) = 2\e) + \f) contains only the two states 
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without arrows (defects) |e) and |/). The various tran- 
sition rates can be obtained from the matrix elements of 
H. 

In Appendix B we are going to use a 2 L dimensional 
representation of the L— 1 generators ej and of the Hamil- 
tonian (IA2|) . In this representation, the Hamiltonian de- 
scribes a spin 1/2 quantum chain. Where can we find 
the eigenvalues of the left ideal (their number is given by 
(|A8p ). among the 2 L eigenvalues of the quantum chain? 
We are going to give an "almost correct" explanation. 
We take again the case L = 4 as an example. If on each 
of the 4 sites of the chain one takes a spin 1/2 representa- 
tion of sl(2), one finds the representation with spin (2 
times), spin 1 (3 times) and spin 2 (one time). If for each 
representation containing 2s + 1 states (s is the spin) one 
takes only the highest weight states, one gets precisely 6 
states (the vertical arrows in Fig. [TOl corresponding to up 
spins) . 

We give now the correspondence between the links- 
defects diagrams and the RSOS configurations considered 
in section 2. For a links-defects diagram with L sites 
(i = 1, 2, . . . , L) take the dual lattice with L + 1 sites (on 
each bond between the sites i and i+1 of the links-defects 
diagram you take the site i on the dual lattice. On the 
dual lattice we have the sites j (j = 0,1,..., L). An 
arrow (defect) on the site i on the links-defects diagram, 
stays unchanged on the dual lattice (it is on the bonds of 
the dual lattice). For the links, one proceeds as follows: 
one takes a site on the dual lattice and a vertical line 
on this site. One counts how many links are cut by the 
vertical line and one takes a vertex with a height h equal 
to the number of intersections. Figs. [2] and [10] illustrate 
the rules. 



In the continuous time limit, the evolution of the sys- 
tem is given by a Hamiltonian H e which corresponds to 
the subspace of highest weight U q (Sl(2) representations 
(q = expi7r/3) [19]. There are ([£/ 2 ]) states in these two 
sectors ([a;] is the integer part of x). If we denote by E r 
(r = 0, 1, . . .) the energy levels in non decreasing order: 
Eq = < Ei < E2 < • • •, the partition function giv- 
ing the finite-size scaling limit of the spectrum of H e is 
defined as follows: 



Z(q) = lim Z L (q) 



hm y q LE n /,v s 



(B4) 



where v s = 3v / 3/2. One can show [2(3] that Z{q) has the 
expression 



(B5) 



Here s is the spin, taking the values s = 0, 1, 2, ... for L 
even and s = |, I , | for L odd, and 



(s(q)=q A °(l-q 2s+1 )fl(l-q n r\ (B6) 



where 



s(2s- 1) 



(B7) 



APPENDIX B: THE FINITE-SIZE SCALING 
LIMIT OF THE HAMILTONIAN 
EIGENSPECTRUM. RESULTS FROM 
CONFORMAL FIELD THEORY. 

We are going to give a brief description of the time 
evolution operator of the stochastic model described in 
Sec. 2. Firstly we consider the spin \ quantum chain 
defined by the Hamiltonian 



(Bl) 



where 



2 to 



v v 



.V5, 
2 



(B2) 



and g x , a v , g z are Pauli matrices. The Hamiltonian (|B1 
commutes with 



"•-IE-*- 



(B3) 



Moreover, for large lattice sizes, the energies are (see 
(dU and (EH)) 



E 



W3 
2L 



(A. + k), 



(B8) 



where k is an integer. 

The Hamiltonian H e has a block diagonal form. The 
states with no defects (L even) and those with one defect 
(L odd) arc in one block. This is the s — (s = h) part 
of (|B5p . The states with defects (L even) and more than 
one defect (L odd) correspond to higher spins. In Sec. 3 
we found that the following values of A s were useful: 



Ai = A 2 = 2 

A 3 = 1, As = — 
3 '23 



(L even) 
(L odd) 



(B9) 



